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ABSTRACT 

We consider how superfluidity of dripped neutrons in tlie crust of a neutron star 
affects the frequencies of the crust's fundamental torsional oscillations. A nonnegligible 
supcrfluid part of dripped neutrons, which do not comove with nuclei, act to reduce the 
enthalpy density and thus enhance the oscillation frequencies. By assuming that the 
quasi-periodic oscillations observed in giant flares of soft gamma repeaters arise from 
the fundamental torsional oscillations and that the mass and radius of the neutron 
star is in the range of 1.4 ^ M/Mq ^ 1.8 and 10 km ^ i? ^ 14 km, we constrain the 
density derivative of the symmetry energy as 100 MeV ^ L < 130 MeV, which is far 
severer than the previous one, L > 50 MeV, derived by ignoring the superfluidity. 
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The structure of neutron stars remains uncertain mainly 
because theoretical extrapolations from empirically known 
properties of atomic nuclei are required to describe the equa- 
tion of state (EOS) of high density matter inside the star. 
Consequently, the theoretical neutron star mass (M) - ra- 
dius (R) relation is still model dependent, whereas the pres- 
ence of a neutron star t hat has mass of abou t 2Mq has 
recently been confirmed (|Demorest et al.l 120101 ) and plays 
a role in ruling out some EOS models. It is generally ac- 
cepted that under an ionic ocean near the star's surface and 
above a fluid core, there exists a crustal region, where nu- 
clei form a bcc Coulomb lattice in a roughly uniform elec- 
tron sea and at a density of about 4 x 10^^ g cm~^, neu- 
trons drip out of the nuclei. The resultant neutron gas is 
a superfluid, as long as the temperature is below the criti- 
cal tempera ture, and presumably relevant to pulsar glitches 
(|Saulslll989l '). In fact, most neutron stars observed are ex- 
pected to be cool enough to have a crust containing a neu- 
tron superfluid. Even at zero temperature, however, a part 
of the neutron gas comoves non-dissipatively with protons 
in th e nuclei via Bragg scattering off the lattice (jChamell 
|2012| ). Crustal torsional oscillations, if occurring in stars, 
would be controlled by the enthalpy density of the con- 
stitu ents that comove with the protons (jvan Horn fc EpsteinI 
Il990l l. In this article, we thus address how neutron su- 
perfluidity affects the oscillation frequencies in a manner 
that depends on the still uncertain EOS of neutron-rich 
nuclear matter. We remark that there are earlier publica- 
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tions that examine the influence of neutron superfluidity on 
the crustal t orsional oscillations for a specific model of the 
cr ust EOS (| Andersson. Glampeda kis. fc Samuelsson' ^2009|; 
Samuelsson fc A ndcrssorJ l2009l : ^Passamonti fc Andcrsso3 
2012h . As we shall see, the effects of neutron superfluidity 



as analyzed here are consistent with the results obtained in 
these publications. 

It is not long ago that oscillations of neutron stars were 
observed in the form of quasi-periodic oscillations (QPOs) 
in giant flares from soft gamma repeaters (SGRs), which 
are considered to be magnetars, i.e., neutron stars with 
surface mag netic fields of order lO""^^ G. This causes 
neutron star asteroseismology to become of practical signif- 
icance in studying the proper ti es of matter in the star (e.g., 
lAndersson fc Kokkotad JlQQahlSotani, Tominaea. fc Maeda 



200]|); ISotani. Kohri. fc Haradal \2004 ): ISotani et al 
(201li)). So far, three giant flares have been detected from 
SGR 0526-66, SGR 1900+14, and SGR 1806-20. The 
frequencies of the QPOs discovered through the timing 
analysis of the a fterglow are in the range from tens Hz 
up to a few kHz (| Watts fc Strohmaverll2006t ). Many theo- 
retical attempts to explain the observed QPO frequencies 
have been done in terms of the torsional oscillations 

in th e crus tal regi on and /or t he magnetic oscillations 

e.g., iLevinI (120061'): iLed (l2007l'l: ISamuelsson fc AnderssonI 



2007h: IS otani. Kokkotas. fc Stergioulaj (200?', '200% 
ISotani. Colaiud a. fc Kokkotas (2008 ) ; Sotani fc Kok kotai 
( 2009h '). The important conclusion is that either the tor- 
sional oscillations or magnetic oscillations dominate the 
excited oscillations in magnetiz ed neutron stars, depending 
on the magnetic field strength jColaiuda fc Kokkotasll201ll : 
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Gabler et al1l201ll .[ 2012al ) and configurations (jCabler et alj 
2012bl ). In the absence of the observational informa- 



tion about the magnetic structure of magnetars except 
the inferred field strengths from the obs erved spindown 
l|Kouveliotou et al]|l99 1; iHurlev et al]|l999l ) , we assume for 
simplicity that the QPOs observed in giant flares are asso- 
ciated with the crustal oscillations. Under this assumption, 
one can probe the properties of nuclear matter in the crust 
such as the density depen dence of the symmetry energy 
and neutron superfluidity jSteiner k, Watt3 20091 : ISotanil 
I2OIII : iGearheart et alJ I2OIII : ISotani et aLlbOia T Once the 
density dependence of the symmetry energy is determined, 
one can address how large the reg i on in which pasta nuclei 
(rods, slabs, etc.) llLorenz et al.l 19931: Ovamatsul Il993l ') 
occur would be, if any i Qvamatsu fc lidal 20071 ). 

To quantify the density dependence of the symmetry 
energy, we use the expansion of the bulk energy per nucleon 
near the saturation point of symmetric matter at zero tem- 
perature l|Lattimeij 1 198 ll ) : 
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where n is the nucleon density, a is the neutron excess, wo, 
no, and Ko denote the saturation energy, saturation den- 
sity, and incompressibility of symmetric nuclear matter, re- 
spectively, and So and L are the parameters characterizing 
the symmetry energy coefficient S{n), i.e., 5*0 = S{no) and 
L = 3no{dS/dn)n=no- The parameters wo, no, and So can 
be well deter mined from the empirical masses and radii of 
stable nuclei jOvamatsu fc lida 2003. 1 , while the parameters 
L and Ko are relatively uncertain. In fact, two of us (K.O. 
& K.I.) constructed the model for the EOS of nuclear mat- 
ter in such a way as to reproduce Eq. ([1} in the limit of 
n — ^ no and a — >■ 0. Then, the optimal density distribution 
of stable nuclei was obtained from the EOS model within 
a simplified version of the extended Thomas-Fermi theory. 
Finally, for given y{= —KoSo/SnoL) and Ko, the values of 
Wo, no, and So were obtained by fitting the charge number, 
mass excess, and charge radius calculated from the optimal 
density distribution to the empirical ones. 

To describe matter in the crust, the Thomas-Fermi 
model was generalized by adding dripped neutrons, a neu- 
tralizing background of electrons, and the lattice energy 
within a Wigner-Seitz approximation l|Ovamatsu fc lidal 
This allows one to obtain the equilibrium nuclear 
shape and size as well as t he crust EOS for variou s sets 
of y and K o - Her e, as in lOvamatsu fc lidal (|2007l ) and 
ISotani et all (|2012l ). we adopt the parameter range < 
L < 160 MeV, 180 sj A'o ^ 360 MeV, and y < -200 
MeV fm'^, which can reproduce the mass and radius data 
for stable nuclei and effectively cover even extreme cases 
ijOvamatsu fc Iida !'2003'l (see Table I in lSotani et"al] (|2012l ') 
for the adopted EOS parameters). 

We turn to calculations of the eigenfrequencies of the 
crustal torsional oscillations. We first prepare the equi- 
librium configuration of the non-rotating neutron stars, 
which is the spherically symmetric solution of the Tolman- 
Oppenheimer-Volkoff (TOV) equations. The corresponding 
solution is described in the following metric with the spher- 
ical polar coordinate: 



In constructing a neutron star from its center, one needs 
the EOS for matter in the core as well as the crust EOS. 
Although the core EOS dominates M and R, even the con- 
stituents of the core remain uncertain. We thus concentrate 
on the crustal region without specifying the core EOS. In do- 
ing so, for various sets of M and R and various models for 
the crust EOS, we inte grate the TOV equa tions inward from 
the star's surface, as in llida fc Satol l| 19971 ) and lSotani et all 
( 2012 ). Hereafter, we shall consider typical values of M and 
R, namely, 1.4 ^ M/Mq ^ 1.8 and 10 km J? 14 km. 
Note that such choice of M and R, which is rather arbitrary, 
encapsulates uncertainties of the core EOS. 

In the presence of torsional oscillations, the elastic- 
ity of the crust works as a restoring force. This elastic- 
ity is characterized by the shear modulus ^, which is an- 
other ingredient of the calculations of the oscillation fre- 
quencies. By assuming that spherical nuclei with charge 
Ze form a bcc Coulomb crystal with number density rii, 
one can approximately describe the shear modulus averaged 
over all directions in the limit of zero temperature as jj, — 
0.11Mni{Ze f /a, where a = (3/47r ni)^/^ is the Wigner-Seit z 
radius (Ogat a fc Ichimaru 1990;: IStrohmaver et al.l Il99ll ). 
As shown in ISotani et al.l ( 20121 )^ this modulus, which 
depends on L mainly through the L dependence of Z 
(jOvamatsu fc lidal [2007h . controls the L dependence of the 
oscillation frequencies. Note that pasta nuclei, if present, 
would contribute to the restoring force in a mann er that 
depends on the shapes (jPethick fc PotekhinI 1199^ ). How- 
ever, the values of L that will be constrained from the 
present analysis sugges t that the region where pasta nuclei 
occur is highly limited (fOvamatsu fc lid a 2007). Thus, as in 
iGearheart et al.l ((SoiJ) and Sotani et al.. (2012), we simply 
consider the torsional oscillations to be confined within the 
region of spherical nuclei. 

Calculations of the frequencies of the torsional modes 
on the spherical equilibrium configuration of the crust as 
obtained above can be performed with high accuracy by 
using the relativistic Cowling approximation, namely, ne- 
glecting the metric perturbations. This is due to the incom- 
pressibility of the modes. With the ^-component of the La- 
grangian displacement vector of a matter element given by 
S."^ ~ y{t, r)dePt/ sm9, where Pe is the ^-th order Legendre 
polynomial, the perturbation equation governing the tor- 
sional oscillations can bo derived from the linearized equa- 
tion of motion as (Schumakcr fc Thornc 1983i) 

Mr^y + [(4 + r<i>' - rA')fj.r + fi'r'']y' 

+ [r'//a;'e-^*-(^ + 2)(^-l)^]e2^3^ = 0. (3) 

Here, we assume y — e'"'3^(r), the prime denotes the deriva- 
tive with respect to r, and H is the enthalpy density defined 
as H = e + p with the pressure p and the energy density e. 
We remark that at zero temperature, the baryon chemical 
potential fib can be expressed as /ij, = H/nt with the baryon 
density ni,. The boundary conditions to be imposed in deter- 
mining the frequency uj of the modes are the zero-traction 
condition at the density where spherical nuclei ceases to be 
present in the deepest region of the crust and the zero-torque 
condition at the star's surface (Schumaker fc Thorne|[T983l : 
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2S(i-) 



dt' + e^^^'-^dr^ + r^ide^ + sm^ 



(2) 



ISotani. Kokkotas. fc StergioulasI 120071 '). 

We now consider the effect of neutron superfiuidity. It 
is generally accepted that at a density of about 4 x lO'^^ g 
cm~'^, neutrons start to be dripped out of nuclei and form 
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a superfluid. Even at zero temperature, however, a compo- 
nent that undergoes Bragg scattering from the bcc lattice of 
the nuclei and thus move non-dissipatively with the nuclei is 
still present. According to recent b and calculation s beyond 
the Wigner-Seitz approximation by ICharnell (I2OI2I ). the su- 
perfluid density, which we define as the density of neutrons 
that are not locked to the motion of protons in the nuclei, 
depends sensitively on the baryon density above neutron 
drip. Since the crustal torsional oscillations are transverse, 
only neutrons that are responsible for the superfluid den- 
sity behave in dependently of the associated dis placements 
of the nuclei ()Pethick. Chamel. fc Reddvl l2010l l. However, 
the enthalpy densit y H in Eg. ([3 1 fully contains the super- 
fluid mass density (|lida fc Bavm[l2002l ). We thus subtract 
the superfluid mass density from H and obtain the effective 
enthalpy density H = H x {A — Ns)/A, where A denotes the 
total nucleon number in a Wigner-Seitz cell, and Ns denotes 
the number of neutrons in a Wigner-Seitz cell that do not 
comove with the nucleus. Finally, by substituting H for H 
in Eq. ([3} , one can determine the frequencies of the torsional 
oscillations in a manner that depends on Ns. Hereafter we 
will assume that Ns comes entirely from the gas of dripped 
neutrons. In fact, although it could come partially from neu- 
trons inside a nucleus, the dripped neutron gas dominates A 
in the entire density region except just above neutron drip. 

We show the results for the frequency of the fundamen- 
tal £ = 2 torsional oscillations, 0^2, calculated for various 
values of Ns/Nd with the number A^^ of dripped neutrons in 
a Wigner-Seitz cell. We remark that Nd — Ns corresponds to 
the number of dripped neutrons bound to the nuclei. First, 
for simplicity, we set the ratio Ns/Nd constant throughout 
the density region above neutron drip. At given M, R, and 
Ns/ Nd, we calculate 0^2 for nine sets of the EOS parameters, 
as m ISotani et al] (|2012l '). Since the calculated 0*2 has a neg- 
ligible dependence on Kq, we construct a fitting formula as 
ot2 = Co — ciL + C2L^ , where co, ci, and C2 are positive and 
adjusted in such a way as to well reproduce the calculations 
for the nine cases. The formulas obtained for Ns/Nd = 0, 
0.2, 0.4, 0.6, 0.8, 1, M = I.SMq, and i? = 14 km are plot- 
ted in Fig. [T] From this figure, one can observe that, due to 
neutron superfiuidity, which acts to enhance the shear speed 
defined as Vs = {/i/Hy^^ , 0^2 becomes up to 66 — 82% larger 
for L ^ 76 A MeV and up to 131% larger for L = 146.1 
MeV than that without neutron superfluidity. Addition- 
ally, in this figure we plot the result for the stellar model 
with the value s of Ns/Nd derived from band calculations in 
IChamell (|2012l ) as "n5i/n^" in Table I. Since these values 
are only of order 10-30% at Ub ~ 0.01-0.4no, the increase 
in ot2 is relativity small and close to the result obtained at 
Ns/Nd = 0.2. Such increase in 0^2 is consistent w i th th e 
results of lAndersson. Glampedakis. fc SamuelssorJ (l2009f): 



ISamuelsson fc Andersson ( 20091 ): Passamonti fc AnderssonI 
l|2012h . 

This increase of 0^2 due to neutron superfluidity would 
make the constraint of L severer if the QPO frequencies 
observed in SGRs come from the torsional oscillations. In 
this case, the lowest frequency in the observed QPOs should 
be equal to or larger than the predicted 0^2, because 0^2 
is the lowest frequency among various eigenfrequencies of 
the torsional oscillations. It is also important to note the 
tendenc y that, as M and R become larger, 0^2 becomes 
smaller jSotani et al.|[20l3 ). Consequently, 0^2 obtained for 
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Figure 1. (Color online) Frequencies of the 1 = 2 fundamental 
torsional oscillations of a star having M = I.8M0 and = 14 
km, ot2, plotted as a function of L. The six solid lines from bottom 
to top correspond to the cases of Ns/Nd = 0, 0.2, 0.4, 0.6, 0.8, 
and 1, while the broken line with symbols shows the result from 
the values of Ns/Nd derived by Chamel (2012). The horizontal 
dot-dashed line denotes the lowest QPO frequency observed from 
SCR 1806-20. 
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Figure 2. (Color online) Lmin as a function of Ns/Nd- The hor- 
izontal broken line corresponds to the result from the band cal- 
culations of Ns/Nd by Chamel. 



M = I.SMq and R = 14 km would determine the lower limit 
of the allowed L from the observed QPO frequencies, as long 
as we confine the stellar models within 1.4 ^ M/Mq ^ 1.8 
and 10 km s; i? < 14 km. 

This lower limit of L, hereafter referred to as Lmin, can 
be obtained from the intersections between the horizontal 
dot-dashed line of 18 Hz and the lines of 0*2 in Fig. [T] For 
constant Ns/Nd, the resultant Lmin is shown in Fig. (2] One 
can observe that the value of Lmin, which is 47.6 MeV in the 
absence of neutron superfiuidity (i.e., Ns/Nd = 0), can be 
as large as 125.9 MeV in its presence (i.e., < Ns/Nd ^ 1). 
In addition, we exhibit Lmin = 55.2 MeV, the result from 
the broken line in Fig. [T] This Lmin, wh ich come s from the 
realistic band calculations of Ns/Nd in ICharnel (|2012l ). is 
expected to give a reliable constraint on L. 

Instead of just considering Lmin, we now proceed to 
obtain a more stringent constraint on L by fitting the 
predicted frequencies of fundamental torsional oscillations 
with different values of I to the low-lying QPO frequen- 
cies observed in S GRs. To t his en d, we use the values of 
Ns/Nd derived by IChamell (|2012l ). In the present analy- 
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Figure 3. (Color online) Comparison of the calculated frequen- 
cies of torsional oscillations of a star having M = 1.4Mq and 
_R = 12 km (sohd lines) with the QPO frequencies observed in 
SGR 1806-20 (dot-dashed lines), where we adopt the Chamel data 
for Ns/Nj. The vertical line corresponds to the value of L that is 
consistent with the observations. 

sis, we focus on the observed QPO frequencies lower than 
100 Hz, i.e., 18, 26, 30, and 92.5 Hz in SGR 1806-20 and 
28, 5 4, and 84 Hz in SGR 1900+14 l|Watts fc Strohmavej 
I2OO6I ). In fact, the even higher observed frequencies would 
be easier to explain in terms of multipolar fundamental 
and overtone frequencies 0. Because of the small interval 
between the observed frequencies 26 and 30 Hz in SGR 
1806-20, it is more difficult to explain the QPO frequen- 
cies observed in SGR 1806-20 than th ose in SGR 1900+14 
l|Sotani. Kokkotas. &: Stergioulaj|2007l ). If one identifies the 
lowest frequency in SGR 1806-20 (18 Hz ) as the fundamen- 
tal torsional oscillation with ^ = 3 as in ISotanl (|201ll ). one 
can manage to explain 26, 30, and 92.5 Hz in terms of those 
with £ = 4, 5, and 15. In the case of the typical neutron 
star model with M = 1.4Mq and R = 12 km, we compare 
the predicted frequencies with the observed ones as shown in 
Fig- El One can observe from this figure that the best value of 
L to reproduce the observed frequencies is L = 127.1 MeV, 
where the calculated frequencies, 18.5 Hz (^ = 3), 24.9 Hz 
{£ = 4), 31.0 Hz {£ = 5), and 90.3 Hz (£ = 15), are within 
less than 5% deviations from the observations. 

Let us now extend the analysis to different stellar mod- 
els and to SGR 1900+14. We find that the QPOs observed 
in SGR 1806-20 can be explained in terms of the eigen- 
frequencies of the same £ within similar deviations even for 
different stellar models except for the case with M — 1.4Mq 
and R — 10 km. The obtained best values of L to explain 
the observations are shown in Fig. [3] by solid lines and filled 
symbols. On the other hand, the low-lying QPOs observed 
in SGR 1900+14 can be similarly explained in terms of the 
fundamental torsional oscillations with £ — 4, 8, and 13; the 
obtained best values of L are shown in Fig.|4]by broken lines 
and open symbols. As a result, the allowed region of L where 
the QPO frequencies observed in SGR 1806-20 and in SGR 
1900+14 are reproducible simultaneously lies in the range 
100 MeV ^ L ^ 130 MeV, as long as the oscillating neutron 



t The possibility to explain the higher QPO frequencies with 
the shear oscillations in hadron-q uark mi xed phase in the core of 
neutron stars is also suggested in lSotaiiT Maruvama. fc Tatsumi 
||2012|') . 



Figure 4. (Color online) Values of L at which the calculated 
frequencies of torsional oscillations agree best with the QPO fre- 
quencies observed in SGR 1806-20 (filled symbols with solid lines) 
and in SGR 1900+14 (open symbols with broken lines), where the 
lines are just eye guides. The circles, squares, and diamonds cor- 
respond to the stellar models with R = 10, 12, and 14 km, respec- 
tively. The gray region denotes the allowed values of L that can 
be obtained from both of the QPO observations in SGR 1806-20 
and SGR 1900+14 by assuming that the corresponding neutron 
stars have mass and radius in the range 1.4 ^ M/Mq ^ 1.8 and 
10 km ^ R ^ 14 km. 

stars have mass and radius ranging 1.4 ^ M/Mq ^ 1.8 and 
10 km ^ J? ^ 14 km. It is interesting to compare this con- 
straint with various ex p erimental constrain ts on L (See, e.g., 
iLi. Chen. fc"Kol (|2008h : lTsang et all (|2009| )). which have yet 
to co nverge but seemingly fa vor smaller L. We remark that 
as m ISteiner fc Watts! l|2009l ). one can also explain the three 
QPOs observed in SGR 1900+14 in terms of the £ = 3, 6, 
9 oscillations with comparable or even better accuracy, but 
the optimal values of L, which are typically of order 80 MeV, 
are too low to become consistent with those based on the 
identification of the four QPOs in SGR 1806-20 as the £ = 3, 
4, 5, 15 oscillations. 

In summary, we have examined how neutron superflu- 
idity affects the fundamental torsional oscillations in the 
crustal region of a neutron star for various EOS and stellar 
models. Recent calculations of the superfluid density allow 
us to determine the increase in the frequencies of fundamen- 
tal torsional oscillations, of which the comparison with the 
QPO frequencies observed in SGRs gives a more stringent 
constraint on the parameter L characterizing the density 
dependence of the symmet ry energy than th e earlier work 
ignoring the superfluidity l|Sotani et al.|[201^ ). If the mass 
and radius of the oscillating neutron stars are observation- 
ally determined, one might be able to justify the assumption 
that the QPOs come from the torsional modes and then to 
flx L with reasonable accuracy. We remark that there could 
be another way of identifying the low-lying QPOs, although 
the 26 Hz QPO in_SGR 1806-20 remains to be clearly iden- 
tifled (ISteiner fc Watt3 12009| ). In these identiflcations, the 
18, 30, 92.5 Hz QPOs in SGR 1806-20 correspond to the 
£ = 2, 3, 10 osciUations, while the 28, 54, 84 Hz QPOs in 
SGR 1900+14 correspond to the £ = 3, 6, 9 oscillations. By 
assuming again that the masses and radii of the oscillating 
neutron stars satisfy 1.4 M/Mq ^ 1.8 and 10 km 7? ^ 
14 km, we obtain a range of L of ~ 60 — 80 MeV in which all 
the six QPO frequencies are reproducible. It would be inter- 
esting to note the difference from the L values of ~ 100 — 130 
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MeV obtained above and, if the 26 Hz QPO in SGR 1806-20 
could be successfully identified as a different type of oscilla- 
tion from the torsional ones, to consider which of these two 
allowed ranges of L is more favored. 
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